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Abstract : The fracture surface of a T^Al-based alloy is studied using 
both an atomic force microscope and a standard scanning electron micro- 
scope. Results are shown to be quantitatively comparable. Two fracture 
regimes are observed. It is shown in particular that the roughness index 
characterizing the small lengthscales regime is equal to 0.5. Furthermore, 
the large lengthscales fractal domain is found to spread over nearly six 
^ decades of lengthscales. 
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I. INTRODUCTION 

Since the pioneering work of B. Mandelbrot and coworkers[l], it has been 
^ shown on all sorts of materials (steels [1, 2], aluminium alloys [6], rocks [3], 
C intermetallic compounds [4, 5], ceramics [7]), using various experimental 
q techniques, that fracture surfaces are self-affine and exhibit scaling proper- 
ties on two [4, 5] or three decades [6] of lengthscales. In most cases, the 
roughness index ( is found to lie around the value 0.8, and it was suggested 
that this could well be a universal value, i.e. independant of the material 
and of the fracture mode [6]. As a matter of fact, it is now believed that 
fracture toughness is not correlated to £, but rather to relevant lengthscales 
measured on the fracture surface [9]. 

As far as metallic materials are considered, however, significantly smaller 
exponents are determined through STM (Scanning Tunneling Microscopy) 
experiments, i.e. for lengthscales lying in the nanometer range. Values of 
the roughness index close to 0.6 for fractured tungstene (regular stepped 
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region), or to 0.5 for graphite, are reported by Milmann and coworkers 
[10, 11]. On the other hand, low cycle fatigue experiments on a steel sample 
have led to a value of ( close to 0.6 [12]. More recently, it was shown that a 
new small lengthscales index could indeed be seen with standard scanning 
electron microscopy (SEM) [8], which was associated to a "quasi-static" 
fracture regime. However, this small lengthscales roughness index, lying 
between 0.4 and 0.5, could not be determined very precisely. It is one of 
the scopes of the Atomic Force Microscopy (AFM) experiments to improve 
the precision on this exponent. 

The upsurge of interest in the problem of crack propagation through 
brittle heterogeneous materials, combined with the progress made in statis- 
tical physics in the understanding of line pinning by randomly distributed 
impurities has led very recently to a few interesting models. In fact, it 
was proposed that the fracture surface could be modelled as the trace left 
behind by a line (the crack front) moving through randomly distributed mi- 
crostructural obstacles. The crack front motion is described by a local non 
linear Langevin equation [13] first written by Ertas and Kardar in a very 
different context [14, 15, 16, 17]. This equation describes a large number of 
regimes, depending on the relative values of the prefactors of the non linear 
terms. Although this first model suffers from some weaknesses [25], it sug- 
gests the existence of two fracture regimes. For crack velocities tending to 
zero, i.e. in the vicinity of the so-called "depinning transition" in the line 
trapping problem, the roughness index (±_ perpendicular to the direction of 
propagation of the crack is predicted to be equal to 0.5. In the dynamic 
case, this exponent is predicted to be also 0.5, except in a small region of 
the parameters space (these parameters being essentially the prefactors of 
the non linearities), where (±_ = 0.75. More generally, for a given velocity, 
there is a crossover at some lengthscale from the "quasi-static" to the "dy- 
namic" behaviour ; the cross-over length decreases quite rapidly with the 
crack velocity (velocity to the power (p =-3). Correlatively, it is expected, 
within the framework of this model, that the crossover length decreases with 
increasing stress intensity factor Kj. 

Another interesting model was proposed by S. Roux et al. [18, 19] for 
the fracture of plastic materials: the fracture surface is expected then to 
be a minimum surface [20], the roughness index of which is lying between 
0.4 and 0.5 [21, 22, 23, 24]. Note that this mechanism, although different 
in its physical content than the "line depinning" one described above, is 
also "quasi-static", since it is based upon an equilibrium model. In this 
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Fig. 1. Sketch of the studied samples. 

la- Crack propagation in mode I in a compact tension specimen. 

lb- Sketch of the broken sample: profiles SEMI and SEM2 observed with a standard SEM on one 
of the two fracture surfaces are located respectively close to zones AFMl/1,2 and AFM2, observed 
on the other fracture surface with an AFM. 



case, however, the crossover between the small and the large lengthscales 
regimes should lie around the plastic zone size, and hence, increase with Kj. 
We are faced then with two different models, which both predict a small 
lengthscales regime characterised by a roughness index 0.5, but in one case 
the crossover length is predicted to decrease with the stress intensity factor, 
and in the second case, it increases with that parameter as the plastic zone 
size. 

In this paper, we report two series of experiments performed with a 
SEM and with an AFM : two regions on a fracture surface, corresponding 
to two different stress intensity factors and, very likely, to different crack 
velocities (although those could not be measured) are analysed on the same 
specimen. The results obtained with the SEM are shown to be quantitatively 
comparable to those obtained with an AFM. It is shown that the small 
lengthscales/ "quasi-static" exponent is indeed very close to 0.5, and that the 
crossover length between this regime and the large lengthscales/ "dynamic" 
regime decreases with the distance to the initial notch, thus confirming 
previous results [8]. In the region far from the notch, the scaling domain is 
found to spread over nearly six decades in lengthscales (5nm-lmm). 



II. EXPERIMENTAL 
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A notched Compact Tension specimen (see Fig. la) of the Super-«2 ahoy 
(T^Al-based) is precracked in fatigue. Fracture is achieved through uni- 
axial tension (mode I) with a constant opening rate (0.2 mm/mn). The 
microstructure of our material is mainly constituted of a<2 needles (~ 1/im 
thick and ~ 20/im long) which break in cleavage, within a j3 matrix, the 
plastic behaviour of which was shown to be important as far as the alloy 
fracture toughness is concerned. One of the two surfaces obtained is electro- 
chemically Ni-Pd plated for SEM observations, while the other one is used 
for AFM. 

Two profiles (SEMI and SEM2, see Fig. lb) located respectively just 
behind the fatigue zone and closer to the edge of the specimen, are obtained 
by subsequently cutting and polishing the sample perpendicularly to the 
direction of crack propagation. These profiles are observed with a scanning 
electron microscope Zeiss DSM 960 at various magnifications, ranging from 
x50 to x3000 or xlO 000, with a backscattered electron contrast. Images in 
256 grey levels are registred and the profiles are extracted by image analysis 
(Visilog 4.1.1). The length of the images is 1024 pixels, and adjacent fields 
(overlapping over fifty pixels with each other) have been explored in order 
to build up profiles of 6 000 and 7 000 points. 

Ten profiles are registred in each of the three different regions, AFM1/1 
and AFM1/2, comparable to SEMI, and AFM2, comparable to SEM2 (see 
Fig. lb). The lengths the profiles are 2.5, 20 and 1 /im, in regions AFM1/1, 
AFM1/2 and AFM2 respectively, with 10 000 points for profiles registered 
in zones AFM1/1 and AFM2, and 20 000 points for those in AFM1/2. 

As it can be seen in Fig. lb, there might be a slight difference of lo- 
calisation between SEMI and AFM1/1 or AFM1/2, and between SEM2 
and AFM2, lying within the millimeter range. However, the experimen- 
tal points registered in regions AFM1/1-AFM1/2 and SEMI on one hand, 
and AFM2-SEM2 on the other hand, nicely collapse on the same curve for 
all the statistical analysis performed (see below). Hence, as it is expected 
within the framework of one or the other of the two models quoted previ- 
ously, the stress intensity factor is the relevant parameter, they do not vary 
much within the considered regions. Note that the overlap region of the two 
techniques extends over two decades (see the inserts of Fig. 2 to 4). 

In order to determine the roughness exponent and the cross-over length, 
three methods are used: the 'variable band width' method, the return prob- 
ability and the spectral method (see [26] for more details). In the case of 
the 'variable band width' method, the following quantity is computed : 
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< Zmax(r) > r =< max{2(r / )} ro<r / <ro+r -min{2(r / )} ro<r / <ro+r > ro oc r c (1) 

where r is the width of the window. Zmax(r) is the difference between the 
maximum and the minimum heights z within this window, averaged over 
all possible origins r of the window belonging to the profile. 

On the other hand, ( can be determined from the scaling of the return 
probability or of the power spectrum. In the former method, the probability 
P(r) that the height z goes back to its initial value z(r = 0) is computed. 
We seek for the points r; which correspond to the same height 2, and build 
up the histogram of the differences r = r; — r^, averaging over all accessible 
values of z ; it was shown [27, 28, 26] that the averaged histogram P(r) 
scales in the following way: 

P(r) oc r~ c (2) 
In the last method, the power spectrum S(k) of the profile [27] is computed: 

S(k) oc k~ l ~ 2C - (3) 

where k is the wave vector. 

The results of the 'variable band width' method, being less noisy, are 
analysed first. 

In the case of the AFM1/1 zone, three power-law regimes can be observed 
on the plot of Fig. 2. The roughness exponent which can be determined at 
short lengthscales is close to 1, which is characteristic of a flat surface. This 
regime will be discussed in the following. 

At larger lengthscales, two regimes can be observed, which correspond to 
the "quasi-static" (with (qs — 0.5) and "dynamic" (with ( ~ 0.8) regimes 
already observed [8] on this material. As in ref.[8], the simplest form of 
the crossover function is chosen, i.e. Zmax(r) is fitted with the sum of two 
power laws: 

Zmax(r) = a0 * ((-^-f 5 + (-^-) 0M ) (4) 

<;qs 4q5 

which allows to define the crossover length ^qs- Exponent (qs — 0.5 is 
measured for lengthscales approximately ranging between lnm and 1/im, i.e. 
on roughly three decades. Finally, at lengthscales larger than (qs — 1/im, 
the "universal" exponent ( ~ 0.8 is recovered. 

In the case of AFM1/2 (see Fig. 3), the same behaviour is observed. 
As could be expected, (qs is the same as in the AFM1/1 zone. However, 
the crossover length separating the "flat" and the quasi-static regimes is 
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Fig. 2. Regions SEMI and AFM1/1. Zmax(r) is plotted versus r on a log-log plot. Note that the 

experimental points obtained with the two techniques gently collapse on the same curve (the region 

of overlap of the two techniques extending approximately from 10 nm to l^um). Three power law 

regimes can be observed. The slope 1 at the origin is indicated by a dotted line. 

Two non linear curve fits of the results are proposed: 

- Zmax{r) = a0((^)° 5 + (^)° 78 ) with aO = 0.12 and ( QS = OAfim. 

Zmax(r) = a0(( j^) 5 + ( j^) 84 ) with a0 = 0.49 and £ QS = 1/mi. 

In each case, error bars are estimated from the scattering of experimental results relative to the 

various micrographs or profiles analysed. 

Insert: Region of overlap between AFM (*) and SEM (). 



significantly larger (ranging from 50 nm to 0.1 /im) than in the previous 
case. This discrepancy will be discussed in the following. 

In the case of AFM2 (see Fig. 4), a shrink of the intermediate lengthscales 
regime is observed. If Zmax(r) is fitted as previously (Eq. (4)), ^qs is shown 
to decrease down to 5 nm. However, a slightly better fit is proposed with 
only one power law (see Fig. 4), with exponent ( = 0.78. The difference 
between the two values of ( which both fit the experimental data gives an 
idea of the error bars to be expected on these exponents. For the first time 
to our knowledge, the universal exponent ( ~ 0.8 is observed on five and a 
half decades of lengthscales . 

It can be seen in Fig. 2 and 3 that a fit with exponents (qs = 0.5 
and ( = 0.78 gives also excellent results in zones SEM1/AFM1/1 and 
SEMI /AFM 1/2, leading then to a crossover length ^qs — 0.1/im. 

The return probability and the power spectrum analysis give results 
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Fig 3 
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Fig. 3. Regions SEMI and AFM1/2. Zmax(r) is plotted versus r on a log-log plot. The slope 1 at 
the origin is indicated by a dotted line. 

- Zmax{r) = aO((j^) 5 + 78 ) with aO = 0.11 and ( QS = 0.1/im. 

Zmax(r) = a0((j^) 05 + ( j^) 84 ) with aO = 0.47 and £ QS = 1/im. 

Insert: Region of overlap between AFM (*) and SEM(). 




Fig. 4. Regions SEM2 and AFM2. Zmax(r) is plotted versus r on a log-log plot. The slope 1 at the 

origin is indicated by a dotted line. 

Two fits of the data are proposed: 

- Zmax(r) = aO * r° 78 with aO = 0.08. 

Zmax(r) = a0((j^)° 5 + (|^)° 84 ) with aO = 10 and ( QS = hnm. 

Note that the "dynamic" regime extends over five and a half decades of lengthscales: 5nm-lmm. 
Insert: Region of overlap between AFM (*) and SEM(). 



7 



10 3 



£ 1 ° 
'E 

>, 

5 

S 10" 5 
< 



Q_ 10 



10- 



Fig 5 
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Fig. 5. Return probability for SEM1+AFM1/1, SEM1+AFM1/2 and SEM2+AFM2. Arbitray units 

for P are chosen in order that the three curves can be put on the same plot. 

SEM1+AFM1/1: 

- P(r) = a0(^^ + (M05)0.5) for f < Q.l^m, and P(r) = ^ for r > 0.1/im 

P(r) = a0(^+ (M25)0 5 ) for r < 1/im, and P(r) = ^ for r > 

SEM1+AFM1/2: 

- P(r) = a0(^ + (^)° 5 ) for r < 0.1pm and P(r) = ^tw for r > 0.1pm 

P{r) = a0(^ + (S^) D 5 ) for r < lpm and P(r) = for r > lpm. 

SEM2+AFM2: 



P(r 



P(r) = a0(^ + (^)° 5 ) for r < 50nm and P(r) = for r > 50nm. 
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Fig 6 
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Fig. 6. Power spectrum for SEM1+AFM1/1, SEM1+AFM1/2 and SEM2+AFM2. Arbitray units 
for S(k) are chosen in order that the three curves can be put on the same plot. 
O SEM1+AFM1/1: 



* SEM1+AFM1/2: 

- S(k) = «0((f0 2 ' 56 + 

5 (fc) = aO((|) 268 +(|) 2 ) 

A SEM2+AFM2: 

_ SW -~ S ( 2 k) = aO((±f^) 2 68 + (±f^) 2 ). 

which are perfectly compatible (see Fig. 5 and 6) with both choices (£ = 
0.78 and £ = 0.84). In the case of the power spectrum, the cross-over 
wave vector kq S is related to ^qs by the relation kg 5 =27r /£qs- This leads 
(region AFM1) to kq S ~ 6/im -1 for £qs=1[i>tii, and to kq S ~ 60/im -1 for 
£,QS = 0.1/im, actually corresponding to the cut-offs between a power-law 
regime with either exponent -2.68 ( — 1 — 2( with ( = 0.84, see Eq. (3), Fig. 
6), or exponent -2.56 ( — 1 — 2( with ( = 0.78, see Fig. 6), and a large wave 
vectors power law decrease with exponent —2 = — 1 — 2( with (qs = 0.5. 
Note that for the computation of the power spectra, only those of the SEM 
profiles containing 6 000 and 7 000 points (magnification xlOO) have been 
used, since profiles with 1 024 points lead to very noisy spectra. 



The above results confirm those obtained in [8], where the intermedi- 



S(*) = a0((f0" 8 + (^) 3 ) 

- -s(k) = ao(C-y^ + C-r) 



III. DISCUSSION 
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ate and large lengthscales regimes were interpreted respectively as "quasi- 
static" and "dynamic". They are also in agreement with the experimental 
results obtained by Milman and coworkers [10, 11] and by Mc Anulty et al. 
[12], and with previously quoted theoretical models [13-17,18-24]. However, 
the "moving line" models [13-17] predict the short lengthscales regimes to 
be independant of the crack velocity [29]. This is clear for the first series of 
experiments [8], and also compatible with the results presented here. 

Finally, it has to be noted that the "quasi-static" regime has only been 
observed on metallic materials [30]. Hence, plasticity might be an impor- 
tant factor for the onset of this regime, either because different fracture 
mechanisms are indeed involved within the plastic zone, or because plastic 
dissipation may slow down crack propagation at small lengthscales. 

On the other hand, quantitative experiments for fracture in fatigue, 
where the crack velocity and the load are measured should allow to clarify 
the meaning of this regime and are currently being performed. In particu- 
lar, it is not clear from this first set of experiments that the crack velocity 
is the relevant parameter to be taken into account. Since fracture does not 
occur at constant load, it is not obvious that region SEMI corresponds to 
a smaller stress intensity factor than region SEM2, and our attempts to 
perform measurements of the load and of the crack velocity in that region 
failed, because fracture occurs too quickly (approximately in 200 ms, with 
a sharp variation of the velocity within 10ms). However, it can be noted 
that an upper bound of the maximum crack velocity could be estimated to 
be less than 1% of the speed of sound in the material: hence, no inertial 
effect has to be taken into account. 

It would also be of great interest to relate ^qs to characteristic lengths 
of the microstructure. Note that ^qs obtained for AFM1/1 is of the order 
of the average thickness of the a<2 needles. However, in a previous study 
of the same material [8], a cross-over length of 10/im could be determined 
for fracture in fatigue. The distribution of the stresses felt by the crack 
front during its propagation depending both on the microstructural disorder 
and on the loading conditions, it is expected that ^qs is linked both to the 
microstructure, and to the local stress intensity factor. Further experiments 
on an aluminium alloy will also be performed to help clarifying that point. 

Finally, the "flat" regime which should exist at shortest lengthscales is 
mixed with a non-physical signal linked to the experimental limitations of 
our AFM. For a given scanning rate, an increase of the number of points 
registered along a profile of constant length requires an increase of the acqui- 
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sition frequency. In order to avoid any mechanical resonance of the piezo- 
electric actuator or of the cantilever, the highest frequencies are limited to a 
few tens of kHz, which imposes a cut-off at short distances. This crossover 
is then expected to be higher for AFM1/2 than for AFM1/1, which is effec- 
tively the case. However, the existence of this non-physical "flat" domain 
slightly influences the statistical characteristics of the "quasi-static" regime 
at small lengthscales. As a consequence, AFM1/2 is more "polluted" than 
AFM1/1. Note also that the importance of this artefact may be more or 
less important depending on the statistical methods used. 

In conclusion, it has been shown that the results obtained with AFM 
are quantitatively compatible with those obtained with standard SEM. The 
simultaneous use of both techniques has allowed, for the first time, an obser- 
vation of the "universal" fracture regime, characterised by a roughness index 
( ~ 0.8, over roughly six decades of lengthscales. At small enough length- 
scales, a "quasi-static" regime spreading over three decades has been ob- 
served, allowing a good determination of its roughness exponent (qs — 0.5. 
The cross-over length ^qs between the two regimes decreases in regions 
farther from the initial notch. 
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